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Abstract 

The Liquid Motion Experiment (LME), designed to 
study the effects of liquid motion in rotating tanks, was 
flown on STS 84. LME was essentially a spin table that 
created a realistic nutation motion of scale-model tanks 
containing liquid. Two spherical and two cylindrical 
transparent tanks were tested simultaneously, and three 
sets of such tanks were employed to vary liquid viscosity, 
fill level, and propellant management device (PMD) design. 
All the tanks were approximately 4.5 in. diameter. The 
primary test measurements were the radial and tangential 
torques exerted on the tanks by the liquid. Resonant 
frequencies and damping of the liquid oscillations were 
determined by sine sweep tests. For a given tank shape, the 
resonant frequency depended on fill level. For the 
cylindrical tanks, the resonances had somewhat different 
frequencies for the tangential axis (0.55 to 0.75 times spin 
rate) and the radial axis (0.73 to 0.78 times spin rate), and 
the tangential axis resonance agreed more closely with 
available analytical models. For the spherical tanks, the 
resonant frequencies were between 0.74 to 0.77 times the 
spin rate and were the same for the tangential and radial 
axes. The damping coefficients varied from about 1 to 
3 percent of critical, depending on tank shape, fill level, 
and liquid viscosity. The viscous energy dissipation rates 
of the liquid oscillations were determined from sine dwell 
tests. The LME energy dissipation rates varied from 0.3 to 
0.5 times the estimates obtained from scaling previous 
ground tests and spacecraft flight data. The PMDs 
sometimes enhanced the resonances and energy dissipation 
rates and sometimes decreased them, which points out the 
need to understand better the effects of PMD on liquid 
motion as a function of PMD and tank design. 

Nomenclature 

d tank diameter 


E dissipation rate of kinetic energy 

g a steady acceleration (such as gravity) 

/ mass moment of inertia of the spacecraft about 

spin 

the spin axis 

/ mass moment of inertia about a transverse axis 

through the center of mass 

K proportionality constant of order unity that 

depends on the spacecraft-tank geometry 

m liquid U( I uid maSS 

p(r, 0, z ) liquid pressure field 

Rq distance from the spin axis to the free surface 
T ri g id computed torque for the empty tank and tank 

structure 

T tot a i amplitude of the measured sinusoidal torque 

a angular acceleration 

A X half-power frequency width 

Y viscous damping coefficient 

X nutation frequency 

p liquid viscosity 

<() phase angle 

p liquid density 

a liquid surface tension 

0 cone angle 

Dq spin rate 

Q in spin motor rotation rate 

^wobble w °bble motor rotation rate 

C 0 j sinusoidal angular velocity of spin table 

(0 unsteady angular velocity 

Introduction 

Background 

The Liquid Motion Experiment (LME) is a flight 
investigation of the characteristics of liquid motions in 
spinning, nutating tanks. The data from the flight tests will 
aid in determining the effects of such motions on the 
stability of spinning spacecraft. 


“This paper is declared a work of the U.S. Government and is not 
subject to copyright protection in the United States.” 
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Spacecraft are made to spin for a number of reasons, 
including gyroscopic stiffness, equal distribution of solar 
thermal loads, and positioning of the liquids in a tank over 
the outlet. Just as for any freely-spinning body, a spinning 
spacecraft can have a nutation motion superimposed on its 
steady spin; this nutation is sometimes called “coning.” 
During nutation, the spin axis of a spacecraft rotates 
around its angular momentum vector, which is fixed in 
inertial space in the absence of any external torques or 
forces. Figure 1 illustrates such a motion for an idealized 
spacecraft containing liquid in two off-axis tanks. 



Fig. 1. — Idealized spinning spacecraft. 
The spin axis rotates around the angular 
momentum vector, setting the liquid into 
motion. 


The angle between the spin axis and the angular 
momentum vector is the cone angle 0, which is a measure 
of the magnitude of the nutation. The angular rate at which 
the spin axis rotates around the angular momentum vector 
is called the nutation frequency X. In a body-fixed 
coordinate system, X is proportional to the spin rate Qq 
according to: 


X = Q 


* spin 


*trans J 


( 1 ) 


where I spin is the mass moment of inertia of the spacecraft 
about the spin axis and I (rans is the mass moment of inertia 
about a transverse axis through the center of mass. The 
spacecraft is assumed to be axisymmetric. Since I trans > 
I in is the common mass distribution for a spinning 
spacecraft, Eq. (1) shows that the nutation frequency is 
smaller than the spin rate, X < Q 0 . 


The nutation motion produces an oscillatory motion 
of the liquids in the tanks. This liquid motion dissipates 
kinetic energy at rate £ as a result of the viscous stresses 
caused by the motion. The dissipation leads to an increase 
in the magnitude of the coning motion [1] as shown by 
Eq. (2): 




( 2 ) 


Since I trans > 7^, the rate of change dQldt of 0 is 
positive; that is, 0 increases with time. Equation (2) 
assumes that the liquid is merely a passive dissipator, 
whose motion does not affect the mass distribution (e.g., 
center of mass location). However, the liquid mass fraction 
of many spacecraft approaches 50 percent of the total mass 
[2,3] and Eq. (2) is then a crude approximation. 
Furthermore, when a resonant frequency of the liquid is 
near the nutation frequency the dynamics of the liquid 
motion couples with the nutation, thereby causing the 
cone angle to change rapidly even in the absence of 
dissipation. Even if there is no dynamic coupling, the 
dissipation £ of a resonant liquid motion is extremely 
large, which, according to Eq. (2), will lead to a rapid 
increase of the cone angle. In extreme cases, the cone 
angle may increase so rapidly that the attitude control 
system cannot maintain the attitude of the spacecraft; 
when this occurs, the spacecraft experiences a disastrous 
“flat spin” [1]. 

For these reasons, it is important to quantify the 
characteristics of the liquid motions in a spinning, nutating 
tank. LME was designed to determine liquid resonant 
frequencies, and energy dissipation rate of resonant and 
nonresonant motions as a function of tank shape (cylinders, 
spheres, and propellant management devices), liquid fill 
level, and liquid properties. 

Overview of Liquid Motions in Spinning. Nutating Tanks 

Liquid in a partially-full spinning tank can oscillate in 
two distinctly different modes: free surface waves, and 
inertial (internal) waves. Of these two, the inertial wave 
mode is the more important because the resonant frequency 
of such modes are in the range that can be excited by 
nutation. 


Free surface waves are similar to the sloshing that 
occurs in non-spinning tanks. In a spinning tank, the 
centrifugal acceleration /^(Qq) 2 (R q is the distance from 
the spin axis to the free surface) is analogous to a steady 
acceleration g (such as gravity), so the resonant frequency 
is proportional to Q^{Rld) 112 where d is the tank diameter 
[4] . Since d < £ 0 and the proportionality constant is greater 
than one, the resonant frequency is greater than the spin 
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rate; in fact, unless the fill level is very small, the resonant 
frequency is greater than twice the spin rate [5]. Thus, the 
resonant frequency of free surface waves are always 
higher than the frequency of the nutation motion that 
excites them. Free surface waves cause an oscillation of 
the liquid center of mass, so in general both a torque and 
a force on the tank are produced by the waves. 

Inertial waves are oscillations of the liquid interior 
and can occur even in the absence of a liquid free surface. 
The oscillations are excited by Coriolis accelerations 
induced by nutation of a spinning tank [6]. For a 
symmetrical tank spinning about its ^-symmetry axis, the 
theory shows, for example, that the liquid pressure field 
p(r,0^)e corresponding to an oscillation frequency of 

X is governed by: 




d 2 p 


dr 2 r 2 W 2 


1- 


4Q 


2^ 


Ip 

\dz 2 


= 0 


0 ) 


When the oscillation frequency X > 2ft 0 , this 
differential equation is “elliptic” and the oscillations are 
the free surface sloshing oscillations mentioned above. 
But when X < 20^, the equation is “hyperbolic” and the 
oscillations are inertial waves [5]. No general theory is 
available for the problem under discussion, but approximate 
analytical models indicate that even for this non- 
axisymmetric geometry all the inertial wave resonances 
occur in the frequency range between 0 and 21^, and the 
most prominent resonances have frequencies less than £2 0 
[7,8,9,10]. Since the motion of the free surface is a 
secondary effect of an inertial wave, the liquid center of 
mass does not necessarily oscillate, in which case the 
liquid motion exerts only a torque on the tank. 


transmitted to the Shuttle Payload General Support 
Computer (PGSC), via the LME Experiment Interface 
Unit and the Data Acquisition System. The LME software 
used to control and monitor the flight tests was resident in 
the PGSC. Data from rotating instrumentation (e.g., load 
cells) are transmitted by an optical link from the spin table 
to the stationary hardware housing, thus eliminating the 
need to transmit data through slip rings. A video camera is 
focused on one tank to provide a visual record of the liquid 
motions. 

Spin and Nutation Rates 

The spin rate and nutation frequency of the spin table 
are provided by two electric motors called the “spin” 
motor and the “wobble” motor. The spin motor housing in 
the final design was mounted on a swash-plate-like 
arrangement tilted at 5° to the wobble motor shaft axis, 
and the tilt angle vector of the spin table rotated at the 
wobble motor rotation rate. The unsteady angular velocity 
w and angular acceleration a of a point on the spin table, 
such as one of the LME tanks, is given in a coordinate 
system rotating with the table by: 


© = & wobble® sin (^ + 40 ( 4a > 

a = Ul wobb i e e cos (Xt + <t>) (4b) 

where <|) is a phase angle that depends on the angular 
location of the table point in question. For comparison, the 
angular velocity and acceleration of an off-axis tank in a 
spinning spacecraft are given by: co = A£sin(ta + <|>) and 
a = A, 2 0cos (Xt + <(>). Thus, if Q wobble = 1, the LME spin table 
simulates the actual motion of a spacecraft. The nutation 
frequency of the spin table was given by: 


The approximate analytical model used to predict the 
resonant frequencies for the LME tanks was based on the 
homogeneous vortex model of [ 1 1 ] . The model was exended 
to incorporate viscous effects was used to predict the 
energy dissipation rates [12]. These predictions are 
compared to the LME test results later in this section. 

Description of Liquid Motion Experiment 

LME Hardware Arrangement 

The LME hardware was designed for the double 
adapter plate of two Shuttle middeck lockers. The 
experiment was actually mounted to a double adapter 
plate of two SpaceHab lockers, but this change in 
accommodations did not require any significant alteration 
to the hardware. The flight hardware is basically a spin 
table that can be made to nutate at a prescribed frequency 
at a fixed cone angle of about 5°. The test tanks are 
cylinders and spheres. The data acquired during a test is 


^ ^ spin ® wobble 

instead of X = & wobhUr The differences between Eqs. (4) 
and (5) and the analogous spacecraft relations are not 
critical because the liquids in the LME tanks are still 
subjected to a spinning, nutating motion, although the 
amplitude of the angular acceleration is not quite the same 
as for a spacecraft. Thus, the LME tests are fully capable 
of determining resonant frequencies and energy dissipation 
rates. 

Load Cells and Tank Support Structure 

The primary quantitative data obtained from the LME 
tests were the torques exerted on the test tanks by the 
contained liquid. These torques were measured by 
extremely sensitive load cells. Two load cell structures 
were used for each tank to provide support for the test tank 
and a capability of sensing both radial and tangential 
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torques. The sensing element of each load cell is a set of 
semi-conductor strain gauges mounted to very thin stainless 
steel beams. The entire weight of the tank support ring 
structure, the test tank, and liquid is supported by the 
tensioned springs that connect these beams and the load 
cell holder brackets The plane of the sensing elements 
passes through the center of mass location of the tank and 
ring support structure. When a torque is exerted about the 
radial axis of the tank from the spin table axis to the tank 
center line, the beams of one of the two load cells deflect 
upward, and the beams of the other load cell deflect 
downward. Conversely, when a torque is exerted about the 
tangential axis of the tank (perpendicular to the radial 
axis), one end of each beam of a load cell deflects upward 
and the other end deflects downward, with an identical 
deflection pattern occurring for the beams of the other load 
cell. Thus, by using an appropriate electrical bridge, the 
load cells can measure radial and tangential torques 
simultaneously. (A radial torque also introduces a 
negligible twist to each beam.) For purely axial or radial 
motions of the tank without a change in the tilt of the tank, 
the strains imposed on the strain gauges effectively cancel 
out; thus, the load cells are insensitive to forces. 

The beam-tensioned spring arrangement provides 
adequate axial support of the tank in the weightless 
environment of flight but not in the 1 -g environment of a 
laboratory, in which the tank and support structure “sag” 
significantly. For that reason, the load cells could not be 
calibrated in the laboratory when they supported a tank. 
Instead, each load cell was calibrated by positioning it 
horizontally to eliminate the dead weight acting on the 
sensing elements, and a fixture was used to load the load 
cell statically over a range of torques. The sensitivities of 
a typical load cell as determined from this procedure are: 
0.4 in.-lb/V (4.52xl0 5 dyne-cm/V) for radial torques and 
0.09 in.-lb/V ( 1 .02x 1 0 5 dyne-cm/V) for tangential torques. 
The difference between the two sensitivities is a result of 
the larger moment arm for radial torques compared to 
tangential torques. The smallest torque imposed on the 
load cells was 0.005 in.-lb (5600 dyne-cm). The calibrations 
were repeated after the flight tests. The differences in pre- 
and post-flight calibration were small, so the pre-flight 
calibration was used to reduce the flight data. On-Orbit 
just prior to a flight data test, the torque offset of each load 
cell was determined by rotating the spin table at the test 
spin rate without nutation. The measured torques (caused 
by slight misalignments, etc.) were subtracted 
electronically from the torques measured during a data 
test, so the recorded test torques should have been due only 
to the nutation motion of the tank and liquid. 


Test Tanks 

LME employed two generic tank shapes, a cylinder 
and a sphere. Figures 2 and 3 show the tank designs 




Fig. 3. — LME spherical tank. 


The tanks were cast from a clear polymeric resin. 
Counterweights were attached to the top or bottom of each 
tank to make the geometric center of the tank coincide 
with the center of mass of the tank-structure combination. 

Propellant Management Devices 

In addition to “bare” spherical and cylindrical tanks, 
two cylindrical tanks and two spherical tanks were fitted 
with generic propellant management devices (PMD’s). 
Figure 4 shows the LME PMD’s. 


NASA/TM — 1998-208479 


4 




The cylindrical tank PMD is a circular plate containing 
a 2.23 in. (5.72 cm) diameter hole. It is located in the tank 
so as to bisect the tank into two equal halves. This kind of 
PMD has been found to damp liquid resonances [13]. The 
PMD wasmadeof 1 mm clear polycarbonate sheet material. 

The spherical tank PMD is atypical vane design. It is 
composed of two large circular vanes that extend over the 
entire tank diameter and four small vanes. In a spacecraft, 
the PMD captures a certain quantity of liquid between the 
large vanes during periods of weightlessness and directs it 
to the small vanes over the tank outlet. For LME, the axis 
of the vane system was oriented along a radial axis such 
that the small vanes were located at the outer diameter of 
the tank; this arrangement preserved the symmetry of the 
liquid orientation in the tank. The PMD was made of 1 mm 
clear polycarbonate sheet. 



Sphere PMD 


Fig. 4. — LME propellant management devices. The 
cylinder PMD bisects a cylindrical tank. The axis of 
the vaned sphere PMD is oriented along the radial 
axis of a spherical tank. 

The cylinder PMD bisects a cylindrical tank. The axis 
of the vaned sphere PMD is oriented along the radial axis 
of a spherical tank. 

LME Test Matrix and Test Procedures 

Tank Sets. Liquid Fill Levels, and Liquid Physical 
Properties 

Three sets of four tanks were used in the flight tests. 
Each set consisted of two cylinders and two spheres. The 
first two tank sets were composed of “bare” tanks, and the 
third set was two cylinders and two spheres containing 
PMD’s. For the first lank set, called Set A, the cylinders 
and spheres alternated around the periphery of the spin 
table (i.e., cylinder, sphere, cylinder, sphere). One of each 


kind of tank was pre-filled to the two-thirds full level and 
the other was filled to the one-third level. The liquid was 
de-ionized water having a viscosity of 0.86 cp. One 
cylindrical tank was in the field of view of the video 
camera, so small neutrally-buoyant particles were mixed 
with its liquid to aid in flow visualization. The tanks also 
alternated around the spin table periphery for the second 
tank set, called Set B; for this tank set, the liquid was de- 
ionized water mixed with 55 percent glycerin (by weight) 
to increase the liquid viscosity to 9.5 cp. For the third tank 
set with the PMD’s, called Set C, the two cylindrical tanks 
were adjacent and the two spherical tanks were adjacent. 
These were filled with pure de-ionized water. One 
cylindrical tank was pre-filled to the one-third full level 
and the other was filled to the two-thirds full level. The 
spherical tanks were also pre-filled to the one-third and 
two-thirds full level. The surface tension of the water was 
about 4x1 O' 4 lb/in. (70 dyne/cm), and the surface tension 
of the water-glycerin mixture was slightly less. The density 
of the water was 0.036 lb/in. 3 (1.0 gram/cm 3 ), and that of 
the water-glycerin mixture was 0.041 lb/in. 3 
(1.14 gram/cm 3 ). All liquid property measurements were 
made at 75 °F (23.9 °C). 

Sine Sweep and Dwell Test Characteristics 

Two general types of tests were employed, each 
having a specific objective. For the first type, the spin rate 
was held constant and the nutation frequency was varied 
over a range; these tests were called Sine Sweep Tests. For 
the second type, the spin rate and the nutation frequency 
were held constant for an extended period of time to obtain 
steady state liquid oscillations; these tests were called Sine 
Dwell Tests. 

Sine sweep tests . — The objectives of the sine sweep 
tests were to determine the resonant frequencies of the 
liquid oscillations and the apparent viscous damping of 
the oscillations. At a resonance, the liquid-induced torque 
sensed by the load cells has a local maximum. Thus, the 
resonant frequencies can be determined by examining the 
torque time histories for maxima. The nutation frequency 
was changed at a rate designed to be sufficiently low such 
that the amplitudes of the resonant torques would be near 
their steady-state values; the sweep rate was small initially 
and increased logarithmically in time, so that most of the 
test time was devoted to the range X < & spin where theory 
indicated that prominent resonances occurred [11]. 

The apparent damping g (fraction of critical damping) 
is determined by measuring the distance on the nutation 
frequency axis at which the torque amplitude has decayed 
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to 1/^2 times the resonant peak amplitude. This damping 
coefficient is useful for equivalent mechanical models of 
an inertial wave oscillation, such as a single degree of 
freedom rotating, nutating disk [10]. 

Sine sweep tests for each tank set were conducted for 
three different spin rates, two (nominally 20 rpm and 
14 rpm) of which produced centrifugal accelerations 
sufficiently high to make the effects of surface tension 
negligible. The third spin rate (nominally 4 rpm) was 
sufficiently small such that the Bond number, defined as 
Bo = pRSlfid/o, was of order unity, where p is the liquid 
density and a is the liquid surface tension. For Bond 
numbers of unity or less, the liquid is effectively in a low- 
gravity environment and the interface is highly curved. 
Thus, the 4 rpm spin rate tests simulated a slowly-spinning 
spacecraft and were meant to quantify the influence of 
surface tension on the resonances. 

Sine dwell tests . — The objective of the sine dwell 
tests was to quantify the energy dissipation rates of the 
liquid oscillations under steady state conditions. For these 
tests, the spin rate and the nutation frequency were held 
constant for up 20 min, to ensure that the liquid oscillations 
reached a steady state condition. Two sine dwell tests were 
conducted for each tank set. 

The energy dissipation rate for each tank is computed 
from the measured torque amplitude and the angular 
velocity of the spin table at the tank location by: 

E i = 0-5( Tfotal -Trigid )«, cos <t> f (6) 

Here T total is the amplitude of the measured sinusoidal 
torque, T rigid is the computed torque for the empty tank 
and tank structure, co • is the sinusoidal angular velocity of 
the spin table (i.e. , the product X0) and is the phase angle 
between the torque and the angular velocity sine waves; 
the phase angle <|> is computed from the timing marks 
recorded in the data stream allowing for the rigid body 
torque. The subscript i indicates a radial or tangential 
component of the parameter. The total energy dissipation 
rate is the sum of the radial and tangential dissipation rates. 
When the phase angle is between ±90°, the energy 
dissipation rate is positive according to Eq. (6). This 
means that the liquid extracts energy from the nutation, 
which is eventually transformed into heat by viscous 
stresses (for a spacecraft, the energy is extracted from the 
spin kinetic energy). Conversely, when the phase angle is 
between 90° and 270°, the energy dissipation rate is 
negative, and the liquid contributes energy to the nutation; 
this is unrealistic, since the liquid oscillations would then 
decay rather than continue in a steady state. 


The rigid body torques needed for Eq. (6) are computed 
from the relation: 

Trigid t i ~ hrans^i — spin “ (?) 

where the + sign is used if i corresponds to the radial 
component and the - is used for the tangential component. 
The angular velocity and acceleration are computed from 
the spin motor and wobble motor rotational rates given by 
Eqs. (4) with 0*5° being the cone angle of the spin table. 

Data Analysis 

Data Reduction 

Torque data were acquired for each sine sweep and 
sine dwell test in the form of a time histories of strain 
gauge voltages. The data were sampled 40 times during 
each spin rate period, which is more than adequate to 
resolve the details of the highest frequency sweep rate or 
spin rate used in the tests. 

All the strain gauge voltages were severely 
contaminated by high and low frequency noise from a 
variety of sources including the wobble motor speed 
reducer. Consequently, it was necessary in the post-flight 
data processing to extract the true signal from each strain 
gauge time history by digital filtering. After filtering, the 
strain gauge signals were combined to compute the 
appropriate torques. For sine sweep tests, all data 
components at frequencies that were somewhat higher 
than the highest nutation frequency and somewhat lower 
than the lowest nutation frequency used in the test were 
filtered. For sine dwell tests, for which the frequency of 
the true signal did not vary with time, the filtering closely 
bounded the test nutation frequency. The filtering yielded 
good results except for the 4 rpm spin rate tests; for these 
low spin rate tests, all the low frequency noise could not 
be completely filtered without affecting the true signal. 

Sine Sweep Test Results 

For a sine sweep test, the nutation frequency was varied 
over a prescribed range by increasing the rotation rate of 
the wobble motor exponentially in time from a starting 
value of near zero to a final value of about twice the spin 
motor rotation rate. As an example, Fig. 5 shows the 
wobble motor rotation rate history for Test A0002, 
for which the wobble motor rotation increased from 
0.2 rad/sec (1.9 rpm) to 3.9 rad/sec (37.6 rpm). Since the 
nutation frequency is the difference between the spin 
motor and wobble motor rotation rates, the nutation 
frequency of this test decreased from an initial value of 
18.1 rpm to zero (when the wobble motor spin rate was 
20 rpm after about 900 sec), and then increased to 1 7.6 rpm 
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Fig. 5. — Wobble motor rotation rate for test A0002 
(20 rpm spin rate). Nutation frequency varies from 
18.1 rpm to 0 rpm to -17.6 rpm. 


in the retrograde direction at the end of the test. 
The prograde motion occupied 900 of the 1200 test sec, 
and the retrograde motion occupied 300 sec. Consequently 
the sweep rate through the retrograde motion was probably 
too fast to excite liquid resonances. 

Typical results . — Figure 6 shows a typical sine sweep 
test result: the radial torque amplitude history for Tank 3 
of Test A0101 (2/3 full cylinder, viscosity = 1 cp, spin 
rate = 14. 1 spin rate) as a function of X/Q 0 . In terms of the 
actual test time, the start of the test corresponds to the 
highest value of X/Q 0 (i.e. , the right hand side of the graph) 
and the end of the test to the lowest (negative) value. The 
slight “jaggedness” of this plot is the result of several 
factors, including some residual noise in the data, but the 


primary causes are: ( 1 ) the large number of nutation cycles 
compressed into the graph, (2) the finite data sampling rate 
of the tests, which sometimes missed the peak of a cycle, 
and (3) the plotting software, which draws a straight line 
from one data point to the next without any “smoothing.” 
Figure 6 indicates that there are two liquid resonances, one 
near X/Qg « 0.5 and a second near XA2 0 » 0.76. (The torque 
peak near X/£2 q * 0.85 is the transient associated with the 
start of the test and does not indicate a liquid resonance).The 
normalized resonance at X/£2 q * 0.5 is much reduced 
compared to the one at XA2 q « 0.76. The higher-frequency 
resonance X/Qq « 0.76 is concluded to be the most 
prominent resonance. In fact, the frequency of this 
resonance agrees well with the predicted resonance 
frequency shown previously. 

The torques measured about the tangential axis also 
displayed resonant liquid oscillations. For example. Fig. 7 
shows the tangential torque response for the same Test 
A010 1 . The resonance at X/Qq « 0.76 is apparent in Fig. 7, 
which is the same value as for the radial torque in Fig. 6, 
but the general shape of the torque is significantly different 
from the radial torque. This difference was caused by the 
offset of the center-of-mass of the liquid along the radial 
axis from the center of the tank (as a result of the liquid 
orientation against the outer wall of the tank); this offset 
superimposed a nonresonant rigid-body-like torque on the 
resonant tangential torque that tended to hide the resonant 
peak. Depending on fill level, the resonant peak was, 


0.10 



-0.15 


-0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 

Nutation Freq/Spin Rate 


Fig. 6.— Normalized radial torque for Tank 1 (2/3 full cylinder) in Test A0101. Most prominent resonance 
occurs for X/Qq ~ 0-76. 
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Nutation Freq/Spin Rate 


Fig. 7.— Normalized tangential torque for Tank 1 (2/3 full cylinder) in Test A0101. Torque amplitude depends 
on magnitude of the spin table motion. 

however, sometimes more evident than in Fig. 7, but without nutation conducted before each sweep test was 

generally it was not possible to determine damping values designed to eliminate such offsets in the recorded sweep 

from the tangential torque responses. responses (by compensating for the offset in the software), 

but apparently nutation introduced an additional 
For some fill levels, the resonant frequency determined uncompensated offset, 

from the tangential torques was significantly less than the 

resonant frequency determined from the radial torques. Summary of cylindrical tank results . — Table I 

The homogeneous vortex model is not capable of predicting summarizes the resonant frequency and damping results 

these differences, since it predicts only one inertial wave as determined from both the radial and tangential torque 

resonance for a given fill level [ 1 1 ].The resonant damping data for all the cylindrical tank tests conducted at 20 and 

of the liquid oscillations was computed from the radial 14 rpm (Bo » 1). For some of the tests, the tangential 

torque responses by the half-power method. For example, torque resonances were not distinct enough to allow the 

from Fig. 6, the peak torque amplitude at the resonant resonant frequency or the damping to be ascertained; these 

frequency X = 0.76Q 0 is about 0.09, measured in moment- cases are indicated in the table by “ — 
of-inertia units. Since the tank rigid body torque is 

negligible compared to the liquid-induced torque, the The radial torque resonances for the l/3rd and 2/3 full 

torque amplitude at the half-power points is 0.09/V2 = tanks occurred for nutation frequencies in the range between 

0.064 moment of inertia units. Therefore, the half-power 0.73Q 0 and 0.78QQ, and the tangential torque resonances 

frequency width from Fig.6 is about AX = 0.03£2q. The occurred in the range between 0.55£2 q and 0.75 Qq. These 

viscous damping coefficient g is thus computed as results are compared to the predictions of the homogeneous 

AXJ2X = 0.03f2o/(1.52f2Q) = 0.02, or 2 percent of the vortex model in Fig. 8. For the l/3rd full tank, the model 

critical damping. This level of damping is about the same predictions compare better to the observed tangential 

as the free-surface sloshing damping for a tank of similar torque resonances than to the radial torque resonances, 

size and shape, which is reasonable since most of the although the predictions are good for both sets of resonances 

energy dissipation occurs at the tank walls for either free for the 2/3rd full tank. For both sets of resonances, the 
surface sloshing or inertia waves. ratio AA2 0 at resonance remained nearly constant when ftg 

D.C. offset . — Many of the measured sinusoidal torque was changed from 14 to 20 rpm, which is in agreement 

responses showed a d.c. offset. The steady spin period, with the model. It should be noted that even though the 
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TABLE I.— SUMMARY OF RESONANT FREQUENCY 2^ AND DAMPING G DATA FOR CYLINDRICAL TANKS 



Radial torque 

Tangential 

torque 

Radial torque 

1 

1 


TEST 

Spin rate 
ipm 

HMHj 

BrnB 

SSI 

HjQS 

■SB 

m 

IPM 

BS3 

m 

EQ9 

ESSfl 

m 

COMMENTS 


Tank 1 (p = 1 cp) 

Tank 3 (p = 10 cp) 


A0002 

20.0 

15.4 

1.8 

14.6 

— 

15.2 

20 

14.0 

— 

resonances are not as prominent as for 
14.1 rpm spin rate 

A0101 

14.1 

10.7 

1.7 

10.1 

— 

10.7 

20 

10.7 

— 

second smaller resonance at A. = 7.0 rpm 
(radial and tangential) 


20.0 

15.4 

— 

11.0 

— 

■EM 

— 

— 


resonances amplitudes are small 


14.1 

10.6 

1.4 

8.0 

— 

■w 

1.7 

8.0 

— 

several smaller resonances in radial axis 


Tank 1 (p = 1 cp) 

Tank 2 (p = 1 cp) 


C0I01 14.1 

10.3 and 
10.9 

17 

10.7 

— 

10.4 

24 

10.6 

— 

tanks contain a PMD 



0.00 0.20 0.40 0.60 0.80 1.00 

Tank Fill Fraction 


Fig. 8. — Resonant frequencies for cylindrical 
tanks at high spin rates. 


radial torque responses occurred over a small range of XI 
Qq ratios, these resonances were not an artifact or resonance 
of the physical setup that occurred at some definite wobble 
motor or spin motor speed. For each resonance the actual 
nutation frequencies were quite different than the wobble 
motor and spin motor rotation rates. The video recordings 
confirmed the existence of liquid resonances. During the 
time when the nutation frequency was sweeping through 
a resonance, the amplitude of the vortex-like swirling 
motion of the liquid oscillation increased dramatically and 
then decreased. None of the test results for Bo = 1 (i.e. , the 
“low gravity” tests) showed any resonances. In fact, the 
actual nutation frequencies used in these tests did not 
exceed about 1 = 0.6Q 0 , so the resonances, which are 
predicted to occur at 1 = 0.65f2 0 and higher (as shown by 
Fig. 8) would not have been excited in the tests. The video 
recordings did show, however, that the liquid interface 
was highly curved for Bo = 1, as expected. 


For the bare tanks, the viscous damping coefficient y 
computed from the radial torque responses is of the order 
of 1 to 2 percent of critical. The damping increased with 
fill level and liquid viscosity, as expected. For the 
cylindrical tanks containing a PMD, the damping was 
considerably larger than for “bare” tanks. This result agrees 
with previous findings for a similar PMD in a cylindrical 
tank with hemispherical ends [13] in which resonant 
liquid oscillations found from drop tower tests were damped 
by the PMD. Furthermore, at least for the 2/3 full tank, the 
PMD, which bisected the tank into equal halves, appeared 
to split the radial-axis resonance into two smaller 
resonances with frequencies slightly larger and slightly 
smaller than for the “bare” tank. 

Summary of spherical tank results . — Liquid resonant 
oscillations were not expected to be prominent in the 
spherical tank sweep tests torques, since such oscillations 
can be excited for a spherical tank only through viscous 
stresses at the tank walls. Nonetheless, some resonances 
were observed, especially in the torque response curves 
for the radial axis. The resonances tended to be more 
highly damped than for the cylindrical tank results, and 
thus smaller in amplitude. The observed resonances again 
tended to be clustered near X = 0.75£2q, which in this case 
are 10 to 20 percent lower than the predicted values, as 
shown in Fig. 9. Evidently, the large viscous coupling of 
the liquid to the tank walls (or perhaps small 
nonaxisy mmetric imperfections in the tanks) was sufficient 
to excite resonances. For the spherical tanks, there was 
little or no difference between the radial and tangential 
torque resonant frequencies. The vane-like PMDs of the 
spherical tanks tended to increase the radial axis liquid 
resonance amplitudes slightly but not the tangential axis 
resonances. The increased torque amplitudes were 
accompanied by a tendency for the damping coefficient to 
be decreased. 
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Fig. 9. — Resonant frequencies for spherical 
tanks at high spin rates. 


Just as for the cylindrical tanks, the tests at Bo = 1 did 
not reveal any resonances, again for the same reasons 
discussed for the cylindrical tanks. The liquid interface 
was highly curved, as was expected. Table II summarizes 
all the resonant frequency and damping results for the 
spherical tank sweep tests conducted at 14.1 and 20 rpm. 

Sine Dwell Test Data 

For a sine dwell test, the filtered torque responses 
were all clean sine waves of constant amplitude. For that 
reason, typical torque plots are not shown. The rotation 
rates of the spin motor and the wobble motor were held 
constant during a sine dwell test for a per* od long enough 
for the liquid oscillations to achieve a steady state condition. 

Cylindrical tanks . — Table III summarizes the torque 
and energy dissipation results for all cylindrical tanks, for 
spin rates of £2 q = 20 and 14. 1 rpm. Results for the “low 


gravity” tests £Jq~ 4 rpm are not shown because the phase 
angle of the torque responses with respect to the tank 
angular velocity could not be determined reliably as a 
result of the noise in these data. 

The torque amplitude results are reasonably consistent 
(i.e., they increase with spin rate and nutation rate). The 
computed energy dissipation rates for Tank 1 of Tank Sets 
A, B, and C are also reasonably consistent and in most 
cases positive, especially for the radial axis torques. 
However, the computed energy dissipation rates for Tank 3 
of Tank Sets A and B are generally negative (phase angles 
greater than 90°), which as was discussed above is 
physically unrealistic. Since the energy dissipation rates 
are small, the phase angle between the torque and the 
angular velocity of the tank tended to cluster around 90°, 
so slight errors in determining the phase angle are probably 
responsible for the negative energy dissipation rates in 
many cases. However, in other cases, the phase angle is 
considerably larger than 90°. The load cell calibrations 
and the sign of the torque response of the load cells for 
Tank 3 were checked after the flight and found to be 
correct, so the reason why phase angles for this tank 
tended to be greater than 90° are not apparent. 

The computed negative energy dissipation rates also 
tended to occur more often for torque responses about the 
tangential axis than for the radial axis. This behavior may 
be related to the fact that liquid resonances found in the 
sine sweep tests were more prominent in the torque 
responses for the radial axis than for the tangential axis. 

Although positive and negative energy dissipation 
rates also occurred for tanks containing PMD’ s, the PMD ’ s 
tended to reduce the energy dissipation rate compared to 
the bare cylindrical tanks. This same trend was observed 
for the damping coefficient g computed from the sine 
sweep test results. The energy dissipation data is plotted as 
a function of spin rate and fill level in Fig. 10. 


TABLE II —SUMMARY OF RESONANT FREQUENCY X^ AND DAMPING G DATA FOR SPHERICAL TANKS 



Radial torque 

Tangential 

torque 

Radial torque 

Tangential torque 


TEST 

Spin rate 
rpm 

ESS 

m 

ESI 

m 

EM 

ESI 

m 

KWH 


COMMENTS 


Tank 2 (p = 1 cp) 

Tank 4 (p = 10 cp) 



20.0 

Wfl 

3 2 

— 

— 

15.4 

2.3 

— 

— 

radial resonances only 

A0101 

14.1 

mi 

1 5 

— 

— 

10.7 

2.0 

10.6 

— 

second smaller resonance at X = 7.0 rpm in 
radial axis 


20.0 

15.4 

— 

— 

— 

15.2 

20 

— 

— 

Tank 4 radial resonance is much more 
prominent 

| B0101 

14.1 


1.9 

— 

— 

TiT» 

1.9 

10.6 

— 

Smaller resonances at 7.3 rpm in radial axis 


Tank 3 (p = 1 cp) 

Tank 4 (p = 1 cp) 


| C0101 14.1 

EfZSKSMEMMSi 

■ESI HIM 

tanks contain a PMD | 
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TABLE ID.— SUMMARY OF ENERGY DISSIPATION RATE DATA FOR CYLINDRICAL TANKS 
W = TANK ANGULAR VELOCITY; T = NET TORQUE, <K = PHASE ANGLE BETWEEN T 
AND F; E - ENERGY DISSIPATION RATE 


TEST 

Spin 

rate, 

rpm 

w, 

rad/sec 

MESM 

■ 

fndU’ 

degree 

in. -lb/sec 

^ t ncnthlt 

in.-lb 

ftuvcntial’ 

degree 

Ftu^cntfe]* 

in.-lb/sec 

A0301, Tank 1 



0.00155 

48.3 


0.00903 

86.9 

27X10- 6 

A0301, Tank 3 

EM 


0.00166 

102.2 


0.00213 

521 

75x10’ 

A0401, Tank 1 


0.149 


24.8 

23x10"* 1 

0.00106 

106.9 

EgQ vm 

A0401, Tank 3 


0.149 

Eaaa 

128.2 


0.00628 

133.1 


A0501, Tank 1 

14.1 


0.00136 

45.5 

27xl0' 5 

0.00089 

95.8 


A0501, Tank 3 

14.1 


0.00117 

98.4 


0.00168 

95.7 


A0601, Tank 1 

14.1 


0.001 12 

320 

4.9x10 5 

0.00019 

139.3 

ESBSESM 

A0601, Tank 3 

14.1 

BtH 

0.00095 

153.0 

-43x10 s 

0.00083 

95.1 



20 

0.114 

0.00133 

69.5 


0.00212 

98.1 

~i.7xia 5 


20 

0.114 

0.00256 

159.9 

BB33E3 

0.00254 

100.5 

-26x10 s 

B0401, Tank 1 

IEOH 

0.149 

0.00302 

80.9 


0.00085 

95.1 


B0401, Tank 3 

20 

0.149 

0.00197 

139.8 

KmWum 

0.00101 

115.2 


B0501, Tank 1 

14.1 

0.056 


58.0 

mm aha 

0.00271 

95.5 


B0501, Tank 3 

14.1 

0.056 

Eels 3 

109.3 

-4.9x10* 

0.00260 

103.5 

m wmum 

B0601, Tank 1 

14.1 

0.103 

0.00113 

6.27 

5.8x10’ 

0.00090 

106.4 


B0601, Tank 3 

14.1 

0.103 

0.00046 

111.4 

-8.6x10* 

0.00106 

106.4 

WwSSum 

C0301, Tank 1 



0.00027 

63.5 

6.8x10* 

0.00228 

56.2 


C0301, Tank 2 

mm 


0.00101 

45.6 

4.0x10’ 

0.00275 

16.2 


C0401, Tank 1 

20 

0.149 

0.00002 

118.3 


0.00094 

102.8 


C0401, Tank 2 

20 

0.149 

0.00075 

107.5 


0.00009 

174.5 


C0501, Tank 1 

14.1 

KTHTM 

0.00048 

38.7 


0.00298 

90.5 


C0501, Tank 2 

14.1 


0.00048 

85.5 


0.00197 

85.7 


C0601, Tank 1 

14.1 


0.00071 

9.01 


0.00089 

74.7 

nggFjja 

C0601, Tank 2 

14.1 


0.00042 

165.6 


0.00031 

79.5 




0.30 0.40 0.50 0.60 0.70 

Liquid Fill Fraction 

Fig. 10.— Energy dissipation rates for cylindrical 
tanks. 

Spherical tanks . — Table IV summarizes the results 
for the spherical tanks, for Qq = 20 and 14 rpm; again, the 
low gravity 4 rpm lest results are not shown. 


energy dissipation rate computed from the tangential axis 
torques are negative. Figure 11 shows the computed 
(positive) energy dissipation rates for the spherical tanks, 
for all the higher spin rate tests. The dissipation for the 
bare tanks tended to increase with spin rate and fill level, 
just as for the cylindrical tanks. For the spherical tanks 
with PMDs, however, the dissipation rate tended to decrease 
with an increase in fill level. This is probably the -esult of 
the specific design of the vane PMD; the PMD blocked 
much of the vortex oscillations for the small fill level, 
whereas the PMD was more open with less blockage for 
the higher fill level. The magnitude of the dissipation rate 
for the spherical tanks with PMDs was sometimes larger 
and sometimes smaller than for the comparable bare tank. 

Comparison of Energy Dissipation Rates to Ground and 
Spacecraft Flight Results 

Many ground-test studies of energy dissipation in 
rotating, nutating tanks have been conducted previously at 
very high spin rates to minimize the effects of gravity on 
the liquid orientation. Most of these studies give energy 
dissipation correlations of the form: 


For the spherical tanks, the energy dissipation rates 
computed for the radial torques are nearly all positive. 
However, just as for the cylindrical tanks, some of the 
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TABLE IV — SUMMARY OF ENERGY DISSIPATION RATE DATA FOR SPHERICAL TANKS 
W = TANK ANGULAR VELOCITY; T = NET TORQUE; <D = PHASE ANGLE BETWEEN T 
AND F; E - ENERGY DISSIPATION RATE 


TEST 

Spin rate, 
rpm 

w, 

rad/sec 

in.-lb 

degree 

Eradla]’ 

in.-lb/sec 

Tuu**otl*I’ 

in. -lb 

degree 

in.-lb/sec 

A0301, Tank 2 

20 

0.114 

0.00426 

60.8 


0.00209 



EEDEBH 

20 

0.114 

0.00047 

95.5 


0.00117 



A0401, Tank 2 

20 

0.149 

0.00414 

6.00 

n 

0.00237 

89.8 

5.9xia 7 

A0401, Tank 4 

20 

0.149 

0.00080 

174.7 


0.00121 

85.8 

6.6 x1a 6 

A0501, Tank 2 

14.1 

0.056 

0.00151 

20.9 


0.00101 

89.6 


A0501, Tank 4 

14.1 

0.056 

0.00041 

9.40 

■ iilTill 

0.00081 

712 


A0601, Tank 2 

14.1 

0.103 

0.00287 

15.7 

gftpgjii 

0.00132 

96.9 

-82x10* 

A0601, Tank 4 

14.1 

0.103 

0.00090 

22.0 

KfeHEtl 

0.00090 

93.2 

-Z6xlO* 

B0301, Tank2 

20 

0.114 

0.00246 

25.5 

|H| 

0.00202 

■a 


B0301, Tank 4 

20 

0.114 

0.00061 

63.9 


0.00151 


Bnrw 

B0401, Tank 2 

20 

0.149 

0.00356 

12 2 



90.8 

-26 xlO* 

B0401, Tank 4 

20 

0.149 

0.00086 

50.6 


IumjFEB 

79.3 

21x10’ 

B0501, Tank 2 

14.1 

0.056 

0.00033 

375 

Kmmm 

0.00149 

KEH 


B0501, Tank 4 

14.1 

0.056 

0.00310 

39.1 

KMjfFSB 

0.00105 


EQE&E9I 

B0601, Tank 2 

14.1 

0.103 

0.00198 

5.08 

BUSS! 

0.00132 

882 

zixia 6 

B0601, Tank 4 

14.1 

0.103 

0.00060 

155 

msBEm 

0.00105 

86.9 

29x1a 6 

C0301, Tank 3 

20 

0.114 

0.00458 

735 


0.00103 

131.9 


C0301, Tank 4 

20 

0.114 

0.00118 

75.2 

mwlHum 

0.00089 

119.7 


C0401, Tank 3 


KYFjn 

0.00560 

Wrga 

irnna 

0.00236 

101.7 

-3.6xl0 5 

C0401, Tank 4 

H 


0.00144 

ial 

WSBBfim 

0.00084 

81.1 

9.6x1a 6 

uijlilfluSSSI 

14.1 

0.056 

0.00201 

38.9 

BISEifl 

0.00169 

97.3 

-60x 10* 

EiiMSSI 

14.1 

0.056 

0.00042 

84.2 

mwatSM 

0.00064 

855 

1.4x 10 6 

C0601, Tank 3 

14.1 

0.103 

0.00187 

25.3 

mkssmm 

0.00145 

91.1 

-15x10* 

C0601, Tank 4 

14.1 

0.103 

0.00022 

77.1 

mZGSxfim 

0.00050 

98.6 

-3.9x10* 



Liquid Fill Fraction 

Fig. 11. — Energy dissipation rates for spherical 
tanks. 


where K is proportionality constant of order unity that 
depends on the spacecraft-tank geometry, \x is the liquid 
viscosity, p is the liquid density, d is the tank diameter, and 
0 is the cone angle [14]. The dependency of E on 0 2 
follows from Eq. (2). For comparison to this ground-test 
correlation, a typical LME case is: 2/3 full spherical tank. 
Assuming K - 1 , Eq. (8) predicts that the LME energy 
dissipation rate is E = 0.00014 in.-lb/sec. The LME 
parameters correspond to Test A0401 , Tank 1, for which 


the measured energy dissipation rate, as shown in Table 
IV is about 0.00023 in-lb/sec. The predicted and measured 
values of E are reasonably close considering that the value 
of K in Eq. (8) is not known for the LME geometry and that 
the correlation is based on data that contains a large 
amount of scatter. 

Both Eq. (8) and the LME data plotted in Figs. 1 0 and 
11 agree in the trend that the energy dissipation rate 
increased with an increase in spin rate or fill level and 
decreased with an increase in nutation frequency. 

Nonproprietary data from flight tests of actual spinning 
spacecraft are difficult to find. However, comparisons of 
Eq. (8) to a considerable number of flight tests for 
“anonymous” spacecraft have been made previously [15]. 
It was found that the ratio of the dissipation rate measured 
in flight to that predicted from Eq. (8) varied from about 
0.3 to about 6. The ratios of the various LME test data to 
Eq. (8) vary in the range of 0.3 to 0.5. Consequently, it can 
be concluded that the LME test results compare with and 
predict the available flight test data reasonably well. 

Summary 

The Liquid Motion Experiment has provided valuable 
data on the interaction between spinning liquids and tanks 
in low gravity. Forces which are not measurable on the 
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ground due to the masking effects of gravity have been 
measured on LME. Liquid resonant frequencies and viscous 
damping coefficients were determined by sine sweep 
tests. The observed resonance frequencies for cylindrical 
tanks depended on fill level and were in the range of 0.73 
to 0.78 times the spin rate for resonances about the radial 
axis and in the range of 0.55 to 75 times the spin rate for 
resonances about the tangential axis. The available 
analytical model predicts only one resonance for a given 
fill level, and this prediction agreed rather closely with the 
tangential axis resonance frequencies observed in the 
tests. The resonances for spherical tanks were in the range 
of 0.74 to 77 times the spin rate and did not vary significantly 
between the tangential and radial axes. Energy dissipation 
rates were determined by sine dwell tests. The LME 
energy dissipation rates varied from 0.3 to 0.5 times the 
estimates obtained from previous ground tests and 
spacecraft flight data. LME has both broadened and 
deepened our understanding of how liquid oscillations 
interact with spinning spacecraft, and will lead directly to 
improved spacecraft design. The LME results also point 
out the need for improved analytical models of such liquid 
oscillations. Areas for further work with LME could 
include testing of specific spacecraft designs and tanks 
whose centerline is on the spin axis. 
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radial axis (0.73 to 0.78 times spin rate), and the tangential axis resonance agreed more closely with available analytical models. For 
the spherical tanks, the resonant frequencies were between 0.74 to 0.77 times the spin rate and were the same for the tangential and 
radial axes. The damping coefficients varied from about 1% to 3% of critical, depending on tank shape, fill level, and liquid viscosity. 
The viscous energy dissipation rates of the liquid oscillations were determined from sine dwell tests. The LME energy dissipation 
rates varied from 0.3 to 0.5 times the estimates obtained from scaling previous ground tests and spacecraft flight data. The PMDs 
sometimes enhanced the resonances and energy dissipation rates and sometimes decreased them, which points out the need to 
understand better the effects of PMD on liquid motion as a function of PMD and tank design. 
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